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In this paper we discuss a formulation of relativistic quantum mechanics that uses Euclidean 
Green functions or generating functionals as input. This formalism has a close relation to quantum 
field theory, but as a theory of linear operators on a Hilbert space, it has many of the advantages of 
quantum mechanics. One interesting feature of this approach is that matrix elements of operators 
in normalizable states on the physical Hilbert space can be calculated directly using the Euclidean 
Green functions without performing an analytic continuation. The formalism is summarized in this 
paper. We discuss the motivation, advantages and difficulties in using this formalism. We discuss 
how to compute bound states, scattering cross sections, and finite Poincare transformations without 
using analytic continuation. A toy model is used to demonstrate how matrix elements of e~^ H in 
normalizable states can be used to construct-sharp momentum transition matrix elements. 

PACS numbers: 21.45+v 



I. MOTIVATION 

In this paper we investigate a framework for constructing relativistic quantum mechanical models of few-degree- 
of-freedom systems that are inspired by an underlying quantum field theory. Our interest is physics at the few GeV 
energy scale. Poincare invariance is an important symmetry at this scale since the energies are comparable to the 
mass scale of hadrons. Most models of systems at these energies are motivated by quantum field theory, but their 
connection with quantum mechanical models of a finite number of degrees of freedom is not straightforward. The 
advantage of a quantum mechanical model of a finite number of degrees of freedom is that the theory is linear and 
it can be in-principle solved, even for strongly interacting systems, with mathematically controlled errors. Bound 
systems of particles present no special problems; they are just point-spectrum eigenstates of the mass operator (rest 
energy) . The construction of a unitary multichannel S operator can be performed using the same methods that are 
used in non-relativistic models. 

One approach that has been successfully used to formulate realistic Poincare invariant quantum mechanical models of 
few-hadron systems is Poincare invariant quantum mechanics pQ. In this approach a dynamical unitary representation 
of the Poincare group is constructed on a few-particle Hilbert space. This approach has been successfully applied 
to treat a number of few-hadron or few quark problems. The virtue of this approach is that it is possible construct 
quantum mechanical models with a finite number of degrees of freedom that have a unitary representation of the 
Poincare group, satisfy a spectral condition, and for fixed numbers of particles satisfy cluster separability. These are 
essentially all of the axioms of quantum field theory, except microscopic locality, which cannot be tested experimentally 
and requires an infinite number of degrees of freedom. One of the disadvantages of this approach is that the models 
do not have a straightforward relation to a Lagrangian field theory. This makes it difficult to use a field theory like 
QCD to systematically improve or constrain the models. 

Because of the difficulties discussed in the previous paragraph it is desirable to explore alternate formulations 
of relativistic quantum mechanics that have a more direct relation to Lagrangian field theory while preserving the 
structure of the underlying quantum theory. The alternative that we pursue in this work is to formulate models that 
are motivated by the standard reconstruction of a quantum theory from the field theory. The Euclidean formulation 
of quantum field theory is convenient for this purpose because (1) is has a direct relation to the action or Lagrangian 
through cither formal path integrals or the Dyson expansion and (2) it is possible to use the Euclidean Green functions 
or generating functional to directly construct the physical Hilbert space using the Euclidean reconstruction theorem. In 
this paper we argue that using this framework it is possible to compute all interesting quantum mechanical observables 
without explicit analytic continuation. 

The proposed models are constructed by restricting the number of degrees of freedom. In passing to a few- 
body model all of the axioms of field theory cannot be preserved. One of the attractive features of the Euclidean 
reconstruction theorem is that locality is an independent requirement that can be relaxed without violating the other 
axioms ; 2J. A feature of this approach is that the dynamics is introduced directly by the model Euclidean Green 
functions or generating functional, rather than in a model Hamiltonian or Lagrangian. One of the challenges of this 
approach is to find a robust framework for modeling Euclidean Green functions or generating functionals with the 
required properties. An advantage is that standard Euclidean field theoretic methods can be used to formulate models. 

The construction of the physical Hilbert space of the field theory from a Euclidean generating functional is discussed 
in the next section. The generating functional for a scalar field is used to illustrate the basic construction. The 
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construction of the Poincare Lie algebra is discussed in the third section. The resulting generators are self-adjoint 
operators on the physical Hilbert space. The construction of one-particle states and the computation of finite Poincare 
transformations on these states is discussed in section four. The one-particle states are used to construct scattering 
states in section five. The computation of finite Poincare transformations on the scattering states are also discussed 
in this section. Everything discussed in sections 2-5 uses only Euclidean generating functionals and Euclidean test 
functions. Section six discuss how the results of sections 2-5 are expressed in terms of Euclidean Green functions. This 
is more practical for formulating models. Models of finite number of degrees of freedom are discussed in section seven 
by considering the structure of models of relativistic nucleon-nucleon scattering. Numerical tests of the proposed 
method to compute scattering observables from matrix elements of e~@ H in normalizable states are given in section 
eight using an exactly solvable non-relativistic model. The results suggest that these methods can be used to compute 
cross sections at the few GeV energy scale. Our conclusions and the outlook for this approach are discussed in the 
last section. 



II. QUANTUM MECHANICS - HILBERT SPACE REPRESENTATION 

To illustrate the construction of the model quantum theory we assume that we are given an Euclidean generating 
functional Z[f] associated with a scalar field. For the purpose of illustration we assume that Z[f] has all of the 
properties that are expected of the generating functional of a scalar field theory. These properties include Euclidean 
invariance, reflection positivity, reality, and cluster separability. 

The generating functional relates the Lagrangian of the field theory to the Hilbert space formulation of the theory. 
In Lagrangian field theory the generating functional is formally the functional Fourier transform of an Euclidean 
path-integral measure with action A[(j>] 

f D\6\e~ A ^ +i ' 1 '^ ^ (i) n ^ (i) n 

j L J n y n y 

n times n times 

where / = /(r, x) is a suitable Euclidean test function and S n is the n-point Euclidean Green function, and S£ is 
the connected n-point Euclidean Green function. In this approach the generating functional Z[f] is the dynamical 
input that replaces the Lagrangian or Hamiltonian. A connection with an underlying Lagrangian or path integral is 
not required, however this connection is an important source of phenomenology, which is the primary motivation for 
developing this formalism. 

In what follows Euclidean space-time coordinates are denoted by x = (r, x). Euclidean invariance of the generating 
functional means that 

Z[f E ,a] = Z[f}. (2.2) 

for 

f Eia (x) = f (E-\x- a)) , (2.3) 

where E is an 0(4) rotation and a is a Euclidean space-time displacement. 

The construction of the physical Hilbert space from Euclidean Green functions was given by Osterwalder and 
Schrader 2 3J. A simpler construction in terms of the generating functional was given by Fr6hlich[3]. We use 
Frohlich's approach to illustrate the construction using the generating functional for a "scalar field" . To construct the 
physical Hilbert space Osterwalter and Schrader select a Euclidean time axis and restrict the space of test functions, 
/(t, x), to real- valued Schwartz functions of four Euclidean variables that vanish for negative Euclidean times; 

5 + = {/(r,x)G5|/(r,x)=0 r < 0}. (2.4) 

We call elements of S+ positive-time test functions. In what follows all test functions will be assumed to be positive- 
time test functions unless stated otherwise. 

Osterwalter and Schrader introduce the Euclidean time-reflection operator O defined by 

e/(r,x) = /(-r,x). (2.5) 

The generating functional, Z[f], is reflection positive if for any finite sequence of real positive-time test functions, 
{fi ■ • • fn} € 6>+, the n x n matrix 

Mij := Z[fi - Qfj] > 0; (2.6) 
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is non-negative. This condition is not automatic. It holds for generating functionals for free fields and for some lattice 
truncations of interacting theories. In general it is a requirement on physically acceptable models. In what follows 
we simply assume this condition is satisfied. 

In Frohlich's construction a dense set of normalizable vectors in the physical Hilbert space is represented by complex 
wave-functionals of the form 



B[4>] 



rib 
3 = 1 



(2.7) 



where bi are complex numbers and fj(r, x) are real positive-time test functions. The argument <j> can be thought of as 
an abstract integration variable. This interpretation is motivated by the path integral representation of the generating 
functional. 

The physical scalar product of two such wave-functionals, B[<p], and 



i<t>(.9k) 



k=l 



is given directly in terms of the generating functional by 



(B\C) 



Tib ,Tl c 

E 

3,k 



b*c k Z[g k - Qfj 



Reflection positivity, (2.6 1, is equivalent to the statement 

(B\B) > 0. 



(2.8) 



(2.9) 



(2.10) 



The physical Hilbert space is obtained by identifying wave functions whose difference has zero norm and completed 
by adding convergent sequences of wave functionals. 

The inner product (2.9) is the physical quantum mechanical scalar product, even though the input only involves 



the Euclidean generating functional and positive-time Euclidean test functions! No analytic continuation is used. An 
explicit illustration of this relationship in given in section six for free scalar particles in eq. (6.4 1 and free spin 1/2 
particles in eq. (6.5 1. 



While the computation of the exact generating functional is equivalent to solving the field theory, models of 
generating functionals are easily constructed. For example, consider the representation of the generating functional 
in terms of connected Euclidean Green functions. If we isolate the contribution to the generating functional from the 
connected n-point Green function (see 2.1 1, 



then the full generating functional is the product 



z[f] = l[z c n [f]- 



(2.11) 



(2.12) 



It follows that the matrix that gives the quantum mechanical scalar product, (12. 61), has the decomposition 



M. 



ij ■- 



Z[fi 



II*:./. «./ II 



(2.13) 



which is an infinite Schur product of the matrices M™- . A sufficient condition for positivity of My is that each M^ 
in the Schur product is positive (this is because the Schur product can be expressed as the restriction of the tensor 
product of positive operators to the diagonal subspace). 



Thus, one strategy for constructing models is to use the representation (2.13), to build models starting with the 



free field two-point function. Z%[f] is reflection positive if it is the generating functional for a generalized free field. 
Including a reflection positive Z%[f] in the Schur product (2.13) would give a generating functional for an interacting 



many-body theory, where the dynamical input is a suitable connected Euclidean 4-point function, 



J model 



[f] = Z c 2 [f]Zt[f] 



Positivity of My means that that 



zt[.h 



(2.14) 



(2.15) 



is a positive matrix for any finite sequence of positive-time test functions. The model can be refined by including 
additional factors, Z^[f]. 
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III. RELATIVISTIC INVARIANCE 



Poincare invariance of a quantum theory implies the existence of a unitary representation of the Poincare group on 
the physical Hilbert space [5] . Equivalently there should be a set of ten infinitesimal generators of the Poincare group 
represented by self-adjoint operators satisfying the Poincare commutation relations. 

The relation between the complex Euclidean group and the complex Poincare group is relevant for constructing 
Poincare generators. To understand the connection between the 4-dimensional Euclidean group and the Poincare 
group consider the matrices 

x= ft-zx-iy\ x _ / it - z x - iy\ , g ^ 

\x + iyt + zj I x + iy it + z J \ ■ J 

A simple calculation shows that det(X) = t 2 — x 2 and det(X) = — (r 2 + x 2 ), which are the Minkowski and Euclidean 
invariant (distances) 2 respectively. Both determinants are preserved under the linear transformations 

X -> X 1 = AXB 1 X -> X' = AXB 1 det(A) = det(B) = 1 (3.2) 

where A and B are complex matrices with unit determinant. In general the pair (A, B) defines a complex Lorentz 
or complex 0(4) transformation. Both (A,B) and (—A,—B) correspond to the same linear transformation of the 
coordinates. In general, the transformed coordinates may become complex, but the determinant remains real and 
unchanged. If B = A* then transformation X' = AXB 1 is a real Lorentz transformation; if A, B 6 SU(2) then the 
transformation X' = AXB 1 is a real 0(4) transformation. 

When A,Be SU(2) the transformation X' = AXB 1 is a six parameter subgroup of the complex Lorentz group. 
While this represents a real Euclidean transformation on the Euclidean Hilbert space (without the Euclidean time 
reversal 0), it defines a complex Lorentz transformation on the physical Hilbert space. It is possible to extract the 
ten Poincare generators by considering infinitesimal forms of these complex Lorentz transformations. 

In order to get a unitary representation of the real Poincare group on the physical Hilbert space the generators 
must be self-adjoint. It turns out that the in the physical scalar product breaks Euclidean invariance in just the 
right way to ensure that all of the Poincare generators are self-adjoint. 

We begin by defining real Euclidean transformations on wave functionals by 

(T(E,a)B) [4>] = J^bje^-^ / j)£>a (x) := f 3 (E~\x- a)) . (3.3) 

3=1 

These transformations leave the generating functional invariant, Z[f] = Z[fE. 3 ]- In general they will not preserve the 
positive time constraint. 

Before we use these transformations to construct Poincare generators on the physical Hilbert space, it is useful to 
note that the wave functionals can also be considered as multiplication operators. For example the operator B[cj>] 
acting on the wave functional C[4>] is a wave functional D[<p] defined by 

D[4>] := B[4>]C[4>\ XXA''''"'"""' "'' • ( 3 - 4 ) 

j=l k=l 

This will become important when we formulate scattering asymptotic conditions. Note that in the wave functional 
representation the order of the Minkowski field operators is determined by the order of the Euclidean time support 
of the positive time Schwartz functions. 

Next we consider the real Euclidean transformations, T(E, a), as complex Poincare transformations on the physical 
Hilbert space. It is useful to treat pure rotations, space translations, Euclidean time translations, and Euclidean 
rotations in planes containing the time axis separately. The Euclidean time reversal operator does not commute 
with the last two transformations. 

We define action of the Poincare generators on the wave functionals considered as operators by: 

[P,S[0]] = - i | i [r(/,(O,a)),S][0] | _ o :=- l f;6 J |-e;^^ o -' ) , (3.5) 
[J ■ h,B[d>]] = -^(T(i?(n,e),0)S)M| e=o := -i£ bj ^^^\ (3.6) 
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where i?(n,£) is an ordinary rotation about the n axis by an angle £. For the Hamiltonian we require that f3 > in 
T(J, (/3, a)) to preserve the positive-time support condition: 

[i7,S[0]] = -^(r(/,(/3,O))S)[0] | ^:=-£6^e;;if'<-' ) . (3.7) 

For the boost generators we first restrict the support of the test functions fj in the wave functionals to a cone 
symmetric about the positive Euclidean time axis that makes an angle < x < 7r/2 with the Euclidean time axis 

5 x ,+ :={/e«S + |/(r,x)=0 taiT 1 ^) > *}. (3.8) 

This ensures that the support condition is preserved for sufficiently small rotations. On these wave functionals we 
consider the Euclidean rotation T(i? e (n, £), 0) in the n — r plane through angle £ < ir/2 — \: 

/ M , fi (r,x) :=/,(r',x') f/ S S x ,+, (3.9) 

with 

t' = tcos(£) — xhsin(^) x' & = Xh cos(£) + r sin(£). (3.10) 

The restrictions on the parameters £ and x ensure that initial and final vectors are in the physical Hilbert space. On 
these vectors the rotationless boost generator is defined by 



[K • n, B[<j>]] = ~(T(R e (h, £),Q)B) [0], <=o := - ]T b^e™"*-™^ (3.11) 



where R e (n,£) is the Euclidean space-time rotation (3.10|. Note the absence of the i in the expressions for H and K. 
This is compensated for by the O that appears in the physical scalar product. 

Direct calculations show that the ten operators H, P,J,K satisfy the Poincare commutation relations and are 
formally Hermitian on the physical Hilbert space. Self-adjointness of H, P,J follows because these operators are 
generators of either one parameter unitary groups or a contractive Hermitian semigroup. The contractive nature 
of time Euclidean time evolution is proved using reflection positivity and positivity properties of the generating 
functional [BJ. This also ensures that Hamiltonian satisfies the spectral condition: 

H > 0. (3.12) 

Self-adjointness of K can be established by verifying that it is the generator of a local symmetric semigroup [7] [5] [5]. 

Matrix elements of the Poincare generators in normalizable states can be expressed directly in terms of the generating 
functional 

a ( N b N c 



(B\P\C) = 



da. 



E E b*c k Z[g k j, {0 ^ - Qfj 



(3.14) 



.j=l k=i 



a=0 



(B\n.J\C) = 



E E c *j d k Z [9k,R(A,O,0 - H /.- 
,3=1 fe=l 



(3.15) 



(B\h ■ K\C) (EE c*d k Z[g ktR<iii 

^ > ~i fe=i 



,0,0 



©/,] 



(3.16) 



G 



The formal Hermiticity of the generators defined above can be deduced from these expressions. For example 

(B\rf\C) = (C\H\B)* = 



d_ 
'd/3 



N b N c 

b 3 C *k Z [-fk,I,(f3,Q) + @9j 

. 3=1 fc=l 



£=0 



d_ 

d/3 



N b N c 

zZ Y h 3 C *k Z \-fk + {@9)j,I,(-(3,0)} 
, j=l k=l 



0=0 



N b N c \ 

Y E b AZ[®9jMm ~fk] = (B\H\C) (3.17) 

K j = 1 k=1 J 0=0 

where we have used reality, Z*[f] = Z[—f], Euclidean invariance, and properties of 0. Hermiticity of the rotationless 
boost generators follows using the same argument. The Euclidean time-reversal operator, 0, plays the role of the 
missing factor of i when integrating by parts. 

The commutation relations can be verified by explicit computation, however they also follow as a direct consequence 
of the relation between complex 0(4) and the complex Lorentz group. 

Matrix elements of e~^ H can also be directly computed in terms of the generating functional: 

N b N c 

(B\e^ H \C) b ^kZ[g k j im - 6/,]. (3.18) 

j=i fe=i 

These matrix elements are needed to compute scattering cross sections and only involve elementary quadratures. 
Matrix elements of the mass Casimir operator can be expressed in terms of the Poincare generators 

M 2 = H 2 - P 2 (3.19) 



(B\M 2 \C) : = 



if- 



+ 



d 2 



d/3 2 da 2 



N b N c 

zZ zZ ^kZbkjip,^ - ©/i]| p=a=0 < 

3=1 fc=l 



(3.20) 



Finally we note that the real Euclidean transformations, T{E 1 a) can be expressed in terms of the Poincare generators 
on the physical Hilbert space by 



-(3H+ia.-T> 



T(i?(n,^),0) = e 



iJ-htp 



T(i? e (n,V),0) 



(3.21) 



Thus Euclidean time evolution and rotations in Euclidean space-time planes look like imaginary time evolution and 
Lorentz transformation with imaginary rapidities. 

The operators defined in ( 3.5|3.6|3~ and 3.11) are self-adjoint operators on the physical Hilbert space that satisfy 
the Poincare commutation relations. Formally they can be exponentiated to give a unitary representation of the 
Poincare group on the physical Hilbert space but, as we will see in the next two sections, this exponentiation is never 
needed. 

The expressions for the matrix elements of all of the Poincare generators in normalizable states ( 3.13||3T6 | are 
directly expressed in terms of the Euclidean generating functional and Euclidean test functions. Analytic continuation 
is not used. 



IV. PARTICLES 



Given the dense set of wave functionals of the form (2.7) and the physical scalar product (2.9) the Gram-Schmidt 
method can be formally used to construct a complete orthonormal set of wave functionals B n \(f>\ : 



(B„\B„ 



(4.1) 
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Since the orthonormal wave functionals are complete, normalizable one-particle states are linear combinations of these 
orthonormal wave functionals with square summable coefficients: 



that are eigenstates of the mass square Casimir operator (3.19) of the Poincare group with eigenvalue A 2 in the point 
spectrum: 

J2(B m \M 2 \B n )b n = J2(B m \(H 2 - V 2 )\B n )b n = X 2 5 mn b n . (4.3) 

n n 

These normalizable states are infinitely degenerate because there is an associated wave packet in the particle's 
momentum and spin. For suitable wave packets these normalizable eigenstates can be decomposed into simultaneous 
eigenstates of mass and linear momentum using translations and Fourier transforms: 

(2 -^e-^ a (T(O,a)* A )[0]. (4.4) 

These wave functionals can be given a plane-wave normalization 

(*A,p'|*A,p>=£(p'-p). (4.5) 

The simultaneous eigenstates of mass and linear momentum can be further decomposed into eigenstates of spin, 
and z-component of spin using 

*aj,p,„M = Z) / dR(T(R,0)* XtR - lv )[<l>]D%,(R) (4.6) 

v^-i Js ( U2 ) 

where the integral is over SU(2) and dR is the SU(2) Haar measure. This projection gives the canonical spin. This 
is the spin measured in the rest frame when the particle is transformed to the rest frame by a rotationless Lorentz 
transformation. Different projections can be used to get states of different helicities or light-front spins. The integral 



in equation (4.6) vanishes if there are no states of mass A and spin j. 



The normalization of the states can be chosen so 

'I'a.,' ,.' ,, 'I'a..,,,-.,, = *(p' - p)<5i'A>- ( 4 -7) 

The state \*&\ ,- )P „) is a single particle state if A is in the discrete spectrum of M. 

Since we started from a linear combination of wave functionals, the single-particle state is formally represented by 
a single-particle wave functional 

tfAj'.p'.p'M- (4.8) 

In general it is not trivial to compute finite Poincare transformations in terms of the generators, however if the 
one-particle state is a non-degenerate state (i.e. the theory has no other particles with the same mass and spin) then 
this state necessarily transforms irreducibly with respect to the dynamical unitary representation of the Poincare 
group. It follows that 

(B\U(A,a)\y^ p ^) = 



(B|V'A,p^,Oy £ ^e-^(P> + l P'• a ^, A1 [A- 1 (^)AA c (^)] (4.9) 
where 

( P y = awp) + aW w A (p) = ^T^ 2 (4.10) 

and A c (?)] is a rotationless Lorentz boost that transforms a frame where a particle of mass A is at rest to one where 
it has linear momentum p. 



Normalizable single-particle states have the form 

*AJ,/M = [dp /foP^AJ'.p^'M (4-H) 



(*Aj >ff |*A,i,/> = / ^ E S*(M,P)/(M,P)- (4-12) 

In this formalism there is no distinction between elementary particles and bound states. They describe particles 
because they have discrete mass eigenvalues. 

What is interesting is that it is possible to construct the Hilbert space, Poincare generators, find single-particle states 
and perform finite Poincare transformations on single-particle states using only the Euclidean generating functional, 
and positive-time test functions without performing any analytic continuation. 



V. SCATTERING THEORY 



In a quantum theory scattering states are solutions of the Schrodinger equation that evolve into asymptotically 
separated non-interacting single-particle states or bound states. In quantum field theory there is no free dynamics 
on the physical Hilbert space, so with the exception of one-particle states, there are no states of non-interacting 
particles on the physical Hilbert space. However, because of cluster properties, there are states that look like states 
of asymptotically separated particles. These states evolve like systems of free particles until the particles get close 
enough to interact. 

One natural framework to formulate scattering asymptotic conditions that is applicable in both quantum mechanics 
and quantum field theory is the two-Hilbert space formulation of scattering [10 . In this framework a separate 
many-particle Hilbert space of non-interacting particles is introduced. This space is used to label the states of the 
asymptotically stable particles. There is a mapping from this asymptotic space to the physical Hilbert space that 
adds the correct description of the internal structure of the particles on the physical Hilbert space when the particles 
are asymptotically separated. 

In the asymptotic Hilbert space composite particles are treated like elementary particles with a given mass, spin, 
and momentum distribution. The internal structure of the composite particle is contained in the mapping to the 
physical Hilbert space. In field theories all particles have internal structure due to their self- interactions. In the 
non-relativistic case the mapping from the asymptotic Hilbert space to the physical Hilbert space has the form 

H I (Ai, J<)Pi, M»> (5-1) 

i 

where the product is a symmetrized tensor product of possibly composite particle states with given momentum and 
spin. Normalizable states in the physical Hilbert space are obtained when this mapping is integrated over square 
integrable functions of the momenta and magnetic quantum numbers of each asymptotically stable particle. In this 
way composite particles are treated as elementary particles with mass A, and spin ji in the asymptotic Hilbert space 
while the mapping adds the composite structure of the asymptotically separated bound states in the physical Hilbert 
space. 

In quantum field theory there are two standard approaches to scattering. The most common is the LSZ treatment 
of scattering, which formulates the scattering asymptotic conditions using interpolating fields that create states from 
the vacuum that have the same quantum numbers as single particle states. It has the advantage that the asymptotic 
conditions can be formulated without solving the one-body problem. The price for this advantage is that weak limits 
must be used to calculate scattering matrix elements. The second approach is Haag-Ruelle scattering [TTJfTS] which 
uses non-local fields that create only one-particle states from the vacuum. Haag-Ruelle scattering is the natural 
generalization of standard quantum mechanical scattering in the field theory setting and it has a natural two-Hilbert 
space formulation[13j 14 . Haag-Ruelle scattering states are defined by norm limits, just like in the non-relativistic 
case. Haag-Ruelle scattering is not commonly used in applications because it requires the solution of the one-body 
problem on the physical Hilbert space as input. In this work the one-body solutions discussed in the previous section 
are used to formulate the Haag-Ruelle asymptotic conditions. 

We begin with a heuristic summary of the two-Hilbert space formulation of Haag-Ruelle scattering in Minkowski 
field theory. For simplicity we consider a scalar field theory with a single-particle state of mass A. To construct the 
mapping from the asymptotic Hilbert space to the physical Hilbert space the Fourier transform of the field, <p(p), is 
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multiplied by a smooth function p(p 2 ) of the square of the Minkowski four momentum that is one when p 2 = —A 2 
(the mass of the asymptotic particle) and identically vanishes when —p 2 is in the rest of spectrum of M 2 . 

The product of the Fourier transform of the field and the test function, cf>p(p) := p(p 2 )<p(p), is then Fourier trans- 
formed back to configuration space. The resulting field 4> p (x) is covariant, but no longer local. It has the property 
that when it is applied to the physical vacuum it creates only the single-particle eigenstate of the mass operator with 
mass A. Because of the multiplication by p(p 2 ), <fi P (x) is a well-behaved operator- valued function of time when it is 
smeared over a test function in three space variables. 

The part of 4> p (x) that asymptotically looks like a creation operator is extracted by taking the linear combination 
of 4> and <t> below: 

A(f,t) := -i J 4> p {x) d f(x)dx (5.2) 

where f(x) is a smooth positive-energy solution of the Klein-Gordon equation. Here smooth means that the Fourier 
transform of the t = solution is a smooth function with compact support in the three momentum. Haag and Ruelle 
show that the TV-particle scattering states in the physical Hilbert space exist and are given by the strong limits 

l*±(/i, •••/*)) = lim A(f N ,t)->>A(f u t)\0). (5.3) 

t— ¥±oo 

Next we express the limit (15. 3| in a two Hilbert space framework that can also be used in the wave functional 



representation of the physical Hilbert space. First we write A(f,t) defined in (5.2 1 by expressing the time derivative 
of the field using the commutator with the Hamiltonian and the time derivative of the Klein-Gordon solution by an 
energy factor: 

A(f,t)\0) = 

dxdpe lHt (i[H,<f> p {0,x)] - z^(p)^(0,x)) e -^ t |0)e^^(P)* +l P- x /(p) (5.4) 



(2tt) 3 / 2 

where /(p) is a test function in the three momentum. Integrating over x gives a partial Fourier transform of the field 



so (5.2) becomes 

A(f,t)\0) = e lHt I ([H,^(0,p)]-o; A (p)^(Q,p)) |0)/(p) e -^W t dp. (5.5) 



The time-dependence is in the quantities e lHt and e ^(p)', where the second factor gives the time dependence of the 
positive-energy solution the Klein-Gordan equatio n. T his is expressed as an operator that acts on the wave packet 



/(p) of a free particle of mass A. It follows that (5.2 1 can be interpreted as a mapping from a dense subset of the 



Hilbert space of square integrable functions, /(p), to the physical Hilbert space 

A(/,<)|0)= e ^%e^ nt |/> (5.6) 
where Hq — oj\ (p) is the energy of the asymptotic particle and 

Ar(p) := ([H, 0,(0, p)] - w A (p)0 p (O, p)) . (5.7) 
By repeating this analysis TV-times the products that appear in the Haag-Ruelle formula, 

HA(h,t)---A(f N ,t)\0), (5.8) 

can be expressed as mappings from an TV-particle subspace of the Fock space of non-interacting particles of mass A 
to the physical Hilbert space: 

II^CM) • • • MfN,t)\0) = e* Ht J A N ( Pl , ■ • • , P iv)e- l( ^^ (p * ))t /i(Pi) • • • MpaO^Pi ■ • • d PJV = (5.9) 

iHt I a ( n —, \ ^ — iHrjt 



e int / A N { Pl ,--- ,PN)e- lUot h(pi)---f N (PN)dpi---dp N (5.10) 
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where 



A N (p u ---, PN ) := Y[ {[H, 0,(0, Pi )] - u x (pi)M°> P^) l°>- 



Equation (5.9 1 has the form 



e iHt A N e- iHot \f) 



(5.11) 



(5.12) 



where Hq — ^2uJ\(pi) is the Hamiltonian for N non-interacting particles of mass A. 
In this notation the Haag-Ruelle theorem, (5.3) has the two-Hilbert space form: 



lim |||tt ± (f)> 

t— >±oo 



e lHt A N e- lHot \i) 



by 



0. (5.13) 
Following what is done in standard quantum mechanical multichannel scattering theory, wave operators are defined 



N±\ 



lim e iHt A N e- lHot \f) 



l*±(f)>- 



t— >±oo 

In the field theory case |12j these wave operators satisfy relativistic intertwining relations 

u(A,a)n N ± = n N± (®c/i(A,o)) 



(5.14) 



(5.15) 



that relate the dynamical representation of the Poincare group with the tensor product of n single-particle irreducible 
representations on the n-particle sector of the asymptotic Fock space. 

The scattering states can be expressed using the representation of the physical Hilbert space in terms of wave 
functionals in section 2. The relevant observation is that 



Mo) 



(5.16) 



is a single-particle state of linear momentum p and mass A. The wave functionals (2.7) of section 2 are vectors in 



the physical Hilbert space, even though they are expressed in terms of Euclidean test functions and the Euclidean 
generating functional. The wave functional, 



defined in the previous section, creates a single-particle state of linear momentum p and mass A. 
Thus, if we make the replacements 



,(p) -> * 



in the two-Hilbert space Haag-Ruelle injection operator, An (5.11), then it becomes the wave functional 

w Ai(Pi)*A 4 ,^,pi,/tj) [4>] 



A N {p 1: (n, ■ ■ - pat, fi N ) [0] := JJ ([H,^ Xi ,^ 



(5.17) 



(5.18) 



(5.19) 



where this functional allows for the possibility of scattering of composite particles with higher spin. 
This leads to the following formal expression for S- matrix elements between normalizable states: 



S 



fi 



lim / A t (p 1 ,m)-/;(PM,iiM)e i£ ^ (p *: 



(^Lf(Pi)Mi) • • • ,PM,MA/)|e 1 \A N {p x ,n x , ■ ■ ■ ,p' N ,/j,' N )}x 



e iE ^ W)t /((Pi,Mi) • ■■f' N (P , N,»'N)I[dp i dp' j 



(5.20) 
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where the scalar product is expressed in terms of the Euclidean generating functional. This expression has the form 

S fi = lim (B(t)\e- 2iHt \C(t)). (5.21) 

t— >oo 



-HHt 



While ( 5.21 1 involves operators that are defined in the wave functional representation, the real time evolution operator, 



, is difficult to calculate in this representation. 



Fortunately this quantity can be replaced by a more easily computable quantity using the Kato-Birman invariance 
principle [TS] [TB] [13] [2] which identifies the limits 



n±|f) := lim e im A N e- iHot \f) = lim e l9{H)t A N e' l9{Ho)t \f) 



(5.22) 



for g{x) in a suitable class of admissible functions provided both limits exist. The content of this result is that in the 
large-time limit the surviving terms correspond to situations where both exponents oscillate in phase, which requires 
that both the dynamical and asymptotic energies are the same. Replacing H and Hq by functions of H resp. Hq does 
not change this result provided the function is increasing, with suitable smoothness. A useful choice for g{x) that is 
in the class of admissible functions is 



9{x) 



-fix 



(3>0. 



For this choice the expressions for the wave operator becomes 



OJf) := lim e- me ~ P " A N e me ~ ?Ho 

n— »-±oo 



f) 



where the time parameter t has been replaced by a dimensionless integer n. 

This means that the expression (5.21 ) for the S matrix elements can be replaced by 



S fi = <*+|¥-> = lim / /r(pi,Mi)---A*( PA /,MM) e - me ^ <E ^ <Pl, x 

n— too 



(5.23) 



(5.24) 



(^m(P1'M1' ' ' ' ; PM,MM)|e 2me SH |4/v(pi,//i!-- - .PjViMjv)) 



x e 



fi (Pi , Mi) • • ■ I'm (Pm> Mm ) II dpi dp J ' 



(5.25) 



The virtue of this expression is that for large fixed n, e 2mx can be uniformly approximated by a polynomial for 
x G [0, 1] 



Jlinx 



P n , e {x)\<e Vxe[0,l]. 



(5.26) 



Because the spectrum of e @ H is in the interval [0, 1] and the approximation is uniform the operator satisfies the same 
inequality 



"line 



Pn,e(e 



< € 



(5.27) 



where the norm on the left is the uniform or operator norm and P„, e (x) and e are the polynomial and error that 
appear in equation (5.26). 



Compared to (5.21) equation (5.25) has the form 

S fi = lim Vd m (n)(B(n)|e^ mH |C(n)). 



(5.28) 



where d m (n) are the coefficients of the polynomial in (5.27). This is useful because matrix elements of powers of 
e -m/3H k e t ween wav e functionals B{n)[4>] and C{n)[4>] can be expressed directly in terms of the generating functional 
using ( |3.18[ ). 

While time-dependent methods are not traditionally used in scattering calculations, they have been used successfully 
in non-relativistic few-body calculations [17] . The advantage of the above formalism is that the entire calculation can 
be performed using only Euclidean methods. 
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In order to calculate sharp- momentum transition matrix elements it is necessary to use narrow wave packets. If the 
transition matrix is sufficiently smooth as a function of momentum, it will factor out of the ^-matrix element, allowing 
one to define scattering observables that do not depend on the details of the wave packet. This is an assumption in 
the standard formulation relating time-dependent and time-independent scattering |18| . For sharp initial and final 
wave packets the on-shell transition matrix elements can be approximated by: 

(pi,/xi, ■ ■•Pat,/xjv|T'|pi,Mi 5 ' ■ -P2,Ai2> 



(B\S\C) - (B\C) (5 2g) 



-2iri(B\6i(p f - Pl )\C}- 

After the wave packets are fixed the limit n — > oo in ( |5.25[ ) can be investigated. For a large enough n the term in 
the limit, (5.25), which has the form 

(B(n)\e- 2me ^ H \C(n)), (5.30) 

will be a good approximation to 

(B\S\C) = (B(0)\S\C(0)). (5.31) 

For this value of n, e ~ 2me ^ can then be uniformly approximated by a polynomial in e~^ H which can be evaluated 
using Euclidean methods. 

(S(n)|e- 2me ^ H |C(n)) «^d m (n)(S(n)| e -" l ^|C(n)). (5.32) 

m 

Combining these three approximations gives an approximation to sharp-momentum transition matrix elements using 
matrix elements of e~ n @ H in normalizable states as input. 

Once the scattering states are known their Poincare transformation properties are determined by equation (5.15) 
and the transformation properties (4.9) of the single-particle states. 

It is useful summarize the steps needed to calculate transition-matrix elements. 

1. Solve the one-body problem. These are eigenstates of the mass-square operator with discrete eigenvalues: $[0]. 

2. Use translational and rotational covariance to construct $((A,j),p,/z)[</>]. 

3. Choose a sufficiently narrow set of single asymptotic particle wave packets /(Pi>/4»). The width must be 
sufficiently narrow to factor the transition matrix elements out of the S'-matrix elements. 

4. Construct the two-Hilbert space mappings 

M*) =11 / ([^, < i > ((A,j)Pn^)]-^A(p l )<i>((A,j)p l ,/i. i ))/(p 4 ,M l )rfp l e me "" <E ^ <P,) (5.33) 

i 

5. Pick a large enough n. 

6. Make a polynomial approximation to e 2mx for x £ [0, 1] 

e 2mx « P 2 „ )£ (z). (5.34) 

7. Calculate 

Sfi = (A(f',n)\P 2n , £ (e-P H )\A(f',n)) (5.35) 



Approximate (pi, ^i, • • • p„, jU„|T|pi, Hi, ■ ■ ■ p 2 , /x 2 ) using (5.29|. 
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The result can be expressed directly in terms of the Euclidean generating functional using (3.18|. 

The discussion above assumes that the starting point is a Euclidean generating functional. For models the require- 
ment that all of these approximations converge are model assumptions that restrict properties of the model generating 
functional or Green functions. These are reasonable requirements, since the model Green functions are modeled after 
the field theoretic Green functions, which are expected to have these properties, at least approximately. 

The usual difficulties of realizing the Poincare symmetry are replaced by the requirement of finding reflection 
positive Euclidean invariant Green functions or generating functionals. It is interesting that the scattering matrix 
constructed above will be unitary, even if the Green function is given perturbatively. 

The proof of the existence of the wave operators in the Haag-Ruelle formulation of scattering depends on properties 
of the field theory that may be violated on approximation. For the case of models the existence of channel wave 
operators can be established by a generalization of Cook's method 19J, which gives the following sufficient condition 
for the existence of N-particle wave operators: 

/>±oo 

/ \\({H,A N }-A N H )e- lHot \f)\\dt < oo. (5.36) 
Jo 

For the simple kinds of two-body models discussed in section VII it is easy to see that the generalized Cook condition 



(5.36) translates into regularity properties of the connected four-point function. This is because 

\\([H,A 2 ]-A 2 H Q )e- iH " t \f)f (5.37) 
is linear in the four-point Euclidean Green function and vanishes when the Green function is replaced by the free 



four-point Euclidean Green function. What controls the integral in (5.37) is the difference between the full and free 
four-point Euclidean Green functions, which is the connected Euclidean four-point function. 

Finally we note that even though the calculation of the scattering observables are based on Euclidean quantities, 
the mechanism that leads to the convergence of the wave operators is in-phase oscillations, not an exponential fall-off. 
The scattering states involve strong limits and replace interpolating fields by fields that create single-particle states 
out of the vacuum. As a result, these calculations are not subject to some of the difficulties encountered in scattering 
calculations based on a Euclidean lattice discretization[20 , however it is necessary to be able to accurately compute 
matrix elements of e~@ H . 



VI. GREEN FUNCTION REPRESENTATION 



While the Euclidean generating functionals provide a concise and elegant description of the theory as well as a 
consistent treatment of the few and many-body problems, the direct Green function approach of Osterwalder and 
Schrader may be more appropriate for constructing phenomenological few-body models. 

In the Green function approach the representation of the physical Hilbert space in terms of Euclidean wave func- 
tionals is replaced by sequences of positive-time support functions of the form: 



B[<t>] -> (x|f) 



( f° \ 

/i(xh) 

/2(X 2 1,X 22 ) 



V 



(6.1) 



/ 



with an inner product that is expressed in terms of multipoint Euclidean Green functions. 



|f) :=^2 J dxx--- dx m dy n ■ ■ ■ dyi9m(® x i ' 1 1 Sx m )G„ i+n (xi, 



i j y n 



,yi)/n(yi, • 



(6.2) 



In the Green function representation the support of /„(yi, • • • , y„) is for < yj < y 2 < • • • . Note that the order of 
the support of the Euclidean times is identical to the order of the fields in the corresponding Minkowski Wightman 
function. Reflection positivity is the condition 



(f|f) > 



(6.3) 



and (g|f) is the physical quantum mechanical scalar product. 
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The relation between the Euclidean and Minkowski scalar products is illustrated for the case of a free Euclidean 
two point function, one for spin zero and one for spin 1/2: 



(f|f) = J ' f{x)G 2 {Qx,y)f{y)d 4 xd 4 y = 



(2 



1 f „ „ a pip-i^-y) 
_y d 4 x rfV 4 p/(x)^^/(y) 



pZ _|_ m z 



(27T) 



I f p -'ipo-(xo+yo)+ip-(x-y) 

rf 4 xdV 4 p/(x) , . j-^-, . ,_ /(y) 



(p° + iw m (p))(p° - iw m (p))" 



3„ I.9(P) 12 



2w m (p) 



> 



(6.4) 



where 



For spin 1/2: 



<?(p) := 



(2^) 3 / 2 



d 4 yf(y)e 



-v m (p)y -ip-9 



(f|f) = J f(xh°G 2 (ex,y)f(y)d 4 xd 4 y = 



(6.5) 



(2 



tt) 4 / 



d 4 xd 4 yd 4 p/(x) e^-^-wy 



jp-(0x-y) m ~ P ' 7e 



P^ + 77?/ 



/(y) 



/ 



^(p)^(p)A 



(6.6) 



where 



A+(P) := 



^m(p) +7°7 - P- m l Q 



(6.7) 



We see clearly that in both case the Euclidean integrals with the Euclidean time reversal are identical to the corre- 
sponding Minkowski scalar products. The equality shows that analytic continuation is not required to compute the 
physical scalar product in the Euclidean representation. 

In the Green-function representation the formulas for the Poincare generators are replaced by 



<x|# If) == {0, i/i(xn), + lAr ) /2(x2i,x 22 ), • • • } 



<9 X 21 ^ x 22 



(6. 



(x|P|f) := {0, -^/i(x„), -i ( -4- + t^t ) / 2 (x2i,x 22 ), • • • } 
dxn V9x 2 i 5x 22 ' 



(6.9) 



9 ( d d 

(x|J|f) :={0,-ixn x —5— /i(xn),-i x 2 i x — ^- + x 22 x — — ) / 2 (x 2 i, x 22 ), • • • } 

9xn V 9X21 9X22/ 



(6.10) 



(x|K|f) :={0, xii 



-0x?i 11 9xii 



/i(xu), 
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* 2i i - + ~ x%2 ^ ] /2!XJI • x ' 2 ! ' ' ! • ,(U: " 



For particles with spin these expression are modified as follows 

dxn) ' ( " " 1 <)xi i 



x ?J- ) ( -«u x J- + s ) ( 6 - 12 ) 



where 

£ = iV^Sie^i e^ K$ ) aa , (6.14) 

and 

B = V p S(e^-P, eT f '') m , (6.15) 

Here 5(i?i,i? 2 ) is a finite dimensional representations of 517(2) x SU(2) associated with the type of field. It is 
constructed by expressing the finite dimensional representation of the Lorentz group 5(A) in terms of SL(2, C) 
matrices S(A,A*) and subsequently replacing A and A* by independent unitary matrices, A and B. 

For the case of fermions the Euclidean time reversal operator also includes a factor 7 for each final particle. 

The formulas summarized in this section are discussed in more detail in [21,. 

VII. FEW-BODY MODELS 

A typical application where relativistic few-body methods are used is elastic nucleon-nucleon scattering. This is 
normally treated using the inhomogeneous Bethe-Salpeter equation. We outline the formulation of this problem in 
the Euclidean quantum mechanical representation. 

We consider a model Green function of the form 



95^ 



/ \ 

5 2 (ex n ,x 12 ) 

54(6X21, 6x 2 2, X 2 3,X 2 4) 

V: ■. j 



(7.1) 



For this model we assume that only 5 2 and 54 are non-zero. Furthermore we assume that 5 2 and 54 are related by 
cluster properties: 

5 4 = 5 2 5 2 + 5 C = S + S c . (7.2) 

The Euclidean Bethe Salpeter kernel is defined by 

S4 1 - S^ 1 = -K. (7.3) 

The structure of 5 2 is determined by covariance up to an unknown Lchmann weight. If the weight is a delta function 
in the mass then this is a free field Euclidean Green function. In this case the one-body solutions that are needed 
to formulate the scattering problem are trivial. If the Lehmann weight also includes some continuous spectrum then 
it is necessary to solve a one-body problem to formulate the scattering asymptotic condition. To do this we take an 
orthonormal set of positive-time test functions and use the Gram-Schmidt method to construct an orthonormal set 

(/nl/m) = / rfxdy/:(6x) 7 °5 2 (x - y)/ m (y) = S mn . (7.4) 



Because the invariant Minkowski Green function is defined with a Dirac conjugate field rather than a Hilbert space 
adjoint, the 7 needs to be eliminated from 5 2 to get the continuation to the Wightman function that serves as the 
kernel of the Hilbert space scalar product. This is achieved by including 7 as the spinor part of the 9 operator: 
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In this basis one-body solutions have the form 



|A)=^ Cn |f„) 



where c„ and A are determined solving the eigenvalue problem for discrete A 2 : 

d 2 d 2 

A 2 C„ = E(gZ? — T^)(/n7°©'5'2/m,/,(r,a))r=a=0C m - 
m 

In this case the Euclidean two-point Green function has the form 



1 

(2^ 



5 2 (x-y) := / d 4 pp(m)dm 

1 P + in 



m -P-7e c ip.( x -y) 



where 



and 



The matrix elements have the form 



«7o e = /3 = 7° = ~7o; 7e = 1 % - 
p(m) = p m 5{m - A) + /0 c (m). 

(/,e 7 °5 2 /) 



(27t) 4 J + m z 



pH/(y) 



where 



t A+(p,m) 
ff (2?r) 3 p ( m ) 9 ( p > m ) d P dm 

a; m (p)+7°7-p-m7° 



2w m (p) 

is the positive energy Dirac projector and 

g(p,m) := / Ae- u '»( p ' Xo - <p '7(x). 



(7.5) 
(7.6) 

(7.7) 

(7.8) 
(7.9) 



(7.10) 



(7.11) 



(7.12) 



Single particle eigenstates of mass, linear momentum and spin are constructed from the mass eigenstates (7.6) 



using 



MP, M) = E c » (^3/2 / E f «( r ' ^" lx - a)e- lR ~ lp x d a D^(R)dR 



The Haag-Ruelle operators A(p,p) are 



A(p )M )[> 



( 27r )3/2 y dp ^Q[3 _W >(P)]^(P'/ i )[ r_ ' ' X_S ]|3=a=o 



(7.13) 



(7.14) 



(7.15) 



The scattering asymptotic states of interest are two-body states. iS-matrix elements in normalizable states can be 
computed using the methods discussed in section 5. Equation (5.29) is replaced by: 

(g\s\f) = 
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lim ]T / g(p[^[) 9 (p' 2 ^' 2 )e me - m ^ + ^ 

n— >oo * * / 



x5 4 (x^x^;x 1 ,x 2 )A(p l!All )[x 1 - 2n/3]A(p 2 , M2 )[x 2 - 2n/3K" e ^ A<P ' l)Wl4)) x 

3(P2, A t 2)ff(Pi, Mi) rf Pi rf P2rfpidp 2 dx 2 dx 2 c!x' 1 dx 2 (7.16) 

Here the reflection positivity is limited to requiring that 7j7 2 0S , 4 is non-negative on products of positive-time test 
functions. 

This illustrates how the approximations discussed above can be implemented in a few-body setting. We note that 
even if S\ is calculated perturbatively, the resulting approximate S matrix will be unitary. 

VIII. SCATTERING TEST 

The scattering computations outlined above and in section V are based on the convergence of a sequence of three 
approximations that are performed in a prescribed order. While they should in principle converge for suitable model 
Green functions, that does not imply that the approximations can be sufficiently well-controlled to give converged 
predictions for reactions at the few GeV energy scale of interest. Since to the best of our knowledge this approach to 
scattering, i.e. computing sharp-momentum transition matrix elements using matrix elements of e~ nl3H in normaliz- 
able states as input, has not even been tested in non-relativistic models, we discuss the implementation of this method 
in an exactly solvable non-relativistic model. This has the advantage that all of the approximations can be compared 
to exact results, and the accuracy of the proposed method can be determined precisely. We consider a non-relativistic 
Hamiltonian with a separable potential that has the range and strength of a nucleon-nucleon interaction. The range 
is fixed by the pion mass while the strength is adjusted to bind two nucleons with the deuteron binding energy. 

The the interaction is taken as a Yamaguchi interaction, with Hamiltonian 



H = ky m -\g)\(g\ (k| ff ) = 5 (k) = ^— - . (8.1) 



The transition matrix elements arc 



with 



(k|t|k')=ff(k)r(— +iO+) 5 (k') (8.2) 



m r , 



^+»0 + ) = - . . mn J I (8-3) 



m n l 



where m n = .94GcV, m n = .14GeV, e& = — 2.24MeV and the coupling constant is determined from these parameters 
by 

777 

A= ^(m^ + \/-m n eb) 2 . (8.4) 

To test the approximations we calculate sharp-momentum transition matrix elements using matrix elements of e~@ H 
evaluated between normalizable states using the methods outlined in section five. While the eigenstates of this model 
can be computed exactly, we had to use the spectral expansion of the Hamiltonian to compute matrix elements of 
e~P H . While this is a complicated construction in the non-relativistic case, it is replaced by an elementary quadrature 



(3.18) in the Euclidean case. 

The first approximation is to extract sharp-momentum transition matrix elements using sufficiently narrow wave 
). For this model exact expressions are available for both the S operator and transition 
em we choose Gaussian wave packets in the relative momenta of the form 



packets in equation (5.29 



matrix. For our test prob 



(j>{k) = Ne~ a{k - ka)2 = Ne- {k - ka)2/k ™ (8.5) 



where TV is a normalization constant and fco is the on-shell momentum. We do not choose a particular direction 
because the interaction is pure s-wave, and the on-shell transition matrix elements at a given energy are given by a 
single complex number. 
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Sharp on-shell transition matrix elements computed exactly and approximately from S'-matrix elements using 



(5.29) are compared as a function of the width, k w of the wave packets. The wave packet widths were deter- 
mined by the requirement that the approximate transition matrix elements agree with the exact transition ma- 
trix elements to an accuracy of less that 0.1%. The results are shown in Table 1 as a function of the relative 
momentum. The first column of Table 1 shows the on-shell momentum ko (center of the Gaussian). The sec- 
ond column shows the value of a used to get the error shown in the fourth column. All of the errors in col- 
umn four round up to .1%. The third column lists k w := 1/s/a for each value of ko and the last column is 
the dimensionless ratio k w /ko- The values of a in column 3 are used in all of the calculations in this section. 

Table 1 



ko 


a 


k 


% error 




[GeV] [GeV" 2 ] 


[GeV] 








0.05 


325000 


0.00175412 


1 


0.035 


0.1 


105000 


0.00308607 


1 


0.030 


0.2 


26000 


0.00620174 


1 


0.031 


0.3 


10500 


0.009759 





1 


0.032 


0.4 


5100 


0.0140028 





1 


0.035 


0.5 


3000 


0.0182574 





1 


0.036 


0.6 


2000 


0.0223607 





1 


0.037 


0.7 


1350 


0.0272166 





1 


0.038 


0.8 


1000 


0.0316228 





1 


0.039 


0.9 


750 


0.0365148 





1 


0.040 


1.0 


600 


0.0408248 





1 


0.040 


1.1 


475 


0.0458831 





1 


0.041 


1.2 


400 


0.05 





1 


0.041 


1.3 


330 


0.0550482 





1 


0.042 


1.4 


290 


0.058722 





1 


0.041 


1.5 


250 


0.0632456 





1 


0.042 


1.6 


210 


0.0690066 





1 


0.043 


1.7 


190 


0.0725476 





1 


0.042 


1.8 


170 


0.0766965 





1 


0.042 


1.9 


150 


0.0816497 





1 


0.042 


2.0 


135 


0.0860663 





1 


0.043 



The last column of Table 1 suggests that a .1% error will generally be obtained if this width is less than 3% of the 
on-shell momentum. This property holds over a wide range of momenta in this model. This is a simple transition 
operator so one may anticipate narrower wave packets are needed for more realistic models. This approximation can 
be improved by decreasing the width of the wave packet; it is the largest source of error in the calculations. In any 
realistic calculation it does not have to be better than the experimental resolution. 



The second step is to approximate S'-matrix elements in these Gaussian wave packets using equation (5.25). It is 
important to first pick the wave packets because the n value needed for convergence depends on the width of the wave 
packet. The approximate quantities are 

(<p\S n \<f>) := J dkdk>(k)e- me ^ fe2/m (ke|e- 2me ^ H |k')e- me_ " fc ' 2/m 0(k') (8.6) 

In Tables 2-8 these quantities are computed using the spectral expansion for H. In these calculations the bound-state 
contribution is not included because it vanishes in the large n limit. Tables 2-8 show the real and imaginary parts 
of matrix elements (<j)\{S n - I)\4>) for different values of n for k = 50, 100, 200, 500, 1000, 1500, 2000 MeV. The exact 
value is given at the bottom of each table. Table 9 shows the values of f3, the product ko x (3 and the n-values 
used in our final calculations. Table 9 suggests that /3 should be chosen so ko x j3 is of order unity. Except for 
the fco = 50MeV case, n — 250 or more gives errors for the n-limits that are smaller than the errors made in the 



factorization approximation, (5.29) 



The n dependence of the real and imaginary part of matrix elements of (4>\(S — I)\4>), computed using (5.25), are 
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plotted as a function of n for different values of fco in figures 1-10. Figures 11 and 12 show how fast the neglected 
bound state contribution to the spectral expansion fall off with n for k$ = 1 GeV. 

Table 2: k = 50[MeV], a = 325000[GcV- 2 ] 



n Re (<l>\(S n - IM Im (4>\(S n - l)\<j>) 



50 


-7 


.62976513315350e-l 


-1 


.52406978178214e- 


•1 


100 


-1 


,33113144491104c+0 


-2 


.75546806155677e- 


•1 


150 


-1 


,67324498184421c+0 


-3 


.50392186517949e- 


•1 


200 


-1 


,83391449191883c+0 


-3 


.86136590065981e- 


•1 


250 


-1 


,89273779093641c+0 


-3 


.99389077283171c- 


•1 


300 


-1 


,90951807884485c+0 


-4 


.03228425149703e- 


•1 


350 


-1 


,91324926322936c+0 


-4 


,04093826445785e- 


•1 


400 


-1 


,91389545311571c+0 


-4 


.04245969269804e- 


•1 


450 


-1 


,91398265491004c+0 


-4 


.04266741575048e- 


•1 


500 


-1 


,91399177708580c+0 


-4 


.04268944214980e- 


•1 


550 


-1 


,91399249452470c+0 


-4 


.04269131167755e- 


•1 


600 


-1 


,91399253056843c+0 


-4 


.04269144047622e- 


•1 


650 


-1 


.91399252976074c+0 


-4 


.04269145893613e- 


•1 



ex -1.91399253060872e+0 -4.04269147714400c-l 



Table 3: k = 100 [MeV], a = 105000[GcV- 2 ] 



n 


Re (<f>\(S n -I)\<f>) 


Im (<f>\(S n -I)\ct>) 




50 


-8. 


.73395186664514c-l 


4.95616337213744e- 


•1 


100 


-1 


.34576615227520c+0 


7.59199494502660e- 


•1 


150 


-1 


.49091126760062c+0 


8.39700869213905e- 


•1 


200 


-1 


.51566533604846c+0 


8. 53352070852317c- 


•1 


250 


-1 


.51799902547669e+0 


8.54631681615040e- 


•1 


300 


-1 


.51811943431498c+0 


8.54697554653043e- 


•1 


350 


-1 


.51812278309620c+0 


8.54699376334219e- 


•1 


400 


-1 


.51812288480017c+0 


8.54699423334768e- 


•1 


450 


-1 


.51812290596551c+0 


8.54699435479409e- 


•1 


500 


-1 


.51812290955424c+0 


8.54699438102000e- 


•1 


550 


-1 


.51812290968227c+0 


8.54699438639971c- 


•1 


600 


-1 


.51812288857871c+0 


8.54699427676778c- 


•1 


650 


-1 


.51812275938123e+0 


8.54699356334623c- 


•1 



ex -1.51812291315971e+0 8. 54699438329052c-! 
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Table 4: fc = 200 [MeV], a = 26000[Gev" 2 ] 



n Rc (<f>\(S n - T)\4>) ImAMSn-jM 



50 


-2 


.08408481834932e- 


•1 


4 


,56550768265380e-l 


100 


-3 


.11945696198071e- 


•1 


6 


.85279276148059e-l 


150 


-3 


,38623394392403e- 


•1 


7. 


.44641333032490e-l 


200 


-3 


.42127100784575c- 


•1 


7. 


.52475859236454e-l 


250 


-3 


.42359208499266e- 


•1 


7. 


.52997722504423e-l 


300 


-3 


.42366821259122c- 


1 


7. 


.53015087494957c-l 


350 


-3 


.42366956571344c- 


1 


7. 


,53015382645461c-l 


400 


-3 


.42366963936415c- 


1 


7. 


.53015400508140e-l 


450 


-3 


,42366965122021c- 


1 


7. 


.53015404697748e-l 


500 


-3 


.42366965247355c- 


1 


7. 


.53015405389574e-l 


550 


-3 


42366962667700e- 


1 


7. 


.53015400608940e-l 


600 


-3 


,42366938358008e- 


1 


7. 


,53015348812648c-l 


650 


-3 


.42366884413759c- 


■1 


7 


,53015247619997e-l 



ex -3.42366967477707c-l 7.53015410457076e-l 



Table 5: k = 500[MeV], a = 3000[Gev" 2 ] 



n Re Im (^|(5 n - I)\4>) 

50 -5.93330385580271e-3 9.65738963854834e-2 
100 -7.12681692349637e-3 1.18415504225673e-l 
150 -7.18129565909453e-3 1.19491569949089e-l 
200 -7.18179706798405c-3 1.19502423190014e-l 
250 -7.18179794641838e-3 1.19502444469971e-l 
300 -7.18179794670870e-3 1.19502444477783e-l 
350 -7.18179794671048e-3 1.19502444477784e-l 
400 -7.18179794671081c-3 1.19502444477784e-l 
ex -7.18179797016073c-3 1.19502444795275e-l 



Table 6: k = l[GeV], a = 600[Gc V - 2 ] 

n Re Im (cf>\(S n - I)\<f>) 

50 -1.47024820732811c-4 1.55922557816223e-2 
100 -1.62726188649875e-4 1.79713868930101c-2 
150 -1.62967714125273c-4 1.80219591282450e-2 
200 -1.62968113934903e-4 1.80220978916775e-2 
250 -1.62968113982642e-4 1.80220979403700e-2 
300 -1.62968113982642c-4 1.80220979403721c-2 
350 -1.62968113982753e-4 1.80220979403720e-2 
400 -1.62968113982975c-4 1.80220979403718e-2 



ex -1.62968113982742c-4 1.80220979403858e-2 
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Table 7: k = 1.5[GeV], a = 250[GcV- 2 ] 



n Re {<l>\(S n -I)\<l>) Im (<l>\(S n - I)\</>) 
50 -1.40242356887477e-5 4.66175982621713e-3 
100 -1.54430995201738e-5 5.52235009412764e-3 
150 -1.54679181726403e-5 5.55112445776045e-3 
200 -1.54679281958447C-5 5.55130663695727c-3 
250 -1.54679280390813c-5 5.55130688968718e-3 
300 -1.54679280390813e-5 5.55130688978369e-3 
350 -1.54679280394143e-5 5.55130688978374e-3 
400 -1.54679280388592e-5 5.55130688978377e-3 
ex -1.54679280242191c-5 5.55130688386830e-3 



Table 8: k = 2.0 [GeV], a = 135[GeV" 2 ] 

n Re (<f>\(S n - IM) Im (4>\(S n - I)\4>) 
50 -2.60094316473225e-6 1.94120750171791e-3 
100 -2.82916859895010e-6 2.35553585404449e-3 
150 -2.83171624670953e-6 2.37471383801820e-3 
200 -2.83165946257657e-6 2.37492460997990e-3 
250 -2.83165905312632e-6 2.37492527186858e-3 
300 -2.83165905257121e-6 2.37492527262432e-3 
350 -2.83165905190508e-6 2.37492527262493e-3 
400 -2.83165905234917e-6 2.37492527262540e-3 



ex -2.83165905227843e-6 2.37492527259701c-3 
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Table 9: Parameters 



kO [GeV] /^GeV" 1 ] 


kO x /? 


n 


0.05 


80.0 


4.0 


630 


0.1 


40.0 


4.0 


450 


0.2 


10.0 


2.0 


470 


0.3 


5.0 


1.5 


330 


0.4 


4.0 


1.6 


235 


0.5 


3.0 


1.5 


205 


0.6 


2.5 


1.5 


225 


0.7 


1.6 


1.2 


200 


0.8 


1.4 


1.12 


200 


0.9 


1.05 


.945 


190 


1.0 


1.0 


1.0 


200 


1.1 


0.95 


1.045 


200 


1.2 


0.9 


1.08 


200 


1.3 


0.85 


1.105 


200 


1.4 


0.7 


0.98 


200 


1.5 


0.63 


0.945 


200 


1.6 


0.57 


0.912 


200 


1.7 


0.5 


0.85 


200 


1.8 


0.45 


0.81 


200 


1.9 


0.42 


0.798 


200 


2.0 


0.4 


0.8 


200 



The third approximation is the polynomial approximation to e l2ne . In this application we shift the potential by 
a constant (the binding energy) so the spectrum of e~^ H is strictly between zero and one. Then it is only necessary 
to find a polynomial approximation to e mx for x E [0, 1]. 

The polynomial approximation is made using the Chebyshev expansion: 

1 N 

f(y) « Pn(v) ■= ^c T (y) + £ c k T k (y) (8.7) 

fe=i 



2 ,2k — lit.. , 2k ~ I ir x , . 

c i = "tt — 7 > /(cos(— — --))cos (j— . (8.8) 

3 N+l^ JK K N+12" VJ N + 12' K ' 

k—l 

The Chebyshev expansion is designed for functions on the interval [—1,1]. To make full use of the information in 
this approximation we approximate e~ 2mx = f(2x — 1) by the polynomial approximation to f(y). This gives the 
polynomial 

i , JL 

f(2e-? H - 1) « -CoT o (2e- 0H - 1) + ^ c k T k (2e~^ H - 1) = P N {e~^ H ) (8.9) 

k=l 

where H' = H + e^. It follows that 

T7 W+1 -RH> , n N + 1 

|e 2 --P J v(x)|<2^ T ^ j |||e 2i " e ^ -P JV ( e -^)|||<2^ T ^ I (8.10) 

Chebyshev polynomials are know to be good approximations to the best uniform polynomial approxima- 
tion [22 . Table 10 shows polynomial approximations to e mx for different values of x £ [0,1], and 
n. The errors are between 10~ n and 10~ 13 for the degree of the polynomial only slightly above n. 
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Table 10: Polynomial convergence 



x n deg poly error % 
0.1 200 150 7.939e+00 
0.1 200 200 3.276c+00 
0.1 200 250 1.925c-ll 
0.1 200 300 4.903e-13 
0.1 630 580 3.573e+00 
0.1 630 630 2.069c+00 
0.1 630 680 5.015c-08 
0.1 630 700 7.456c-ll 
0.5 200 150 1.973c-13 
0.5 200 200 1.627e-13 
0.5 200 250 3.266e-13 
0.5 630 580 1.430c-14 
0.5 630 630 3.460c-13 
0.5 630 680 9.330c-13 
0.9 200 150 7.939c+00 
0.9 200 200 3.276c+00 
0.9 200 250 1.950c-ll 
0.9 200 300 9.828e-13 
0.9 630 580 3.573e+00 
0.9 630 630 2.069c+00 
0.9 630 680 5.015c-08 
0.9 630 700 7.230c-ll 



Table 10. suggests that a polynomial with degree 10-20% larger than n is needed for convergence to one part in 10 12 

(io- 10 %). 

All three approximations arc combined to get the approximations to the on-shell transition operator for inci- 
dent momenta between 50 MeV and 2 GeV. These results are displayed in table 11. The errors are all bet- 
ter than .1%. The only significant source of error is the approximate factorization of the sharp momentum 
matrix element from the S matrix elements. This can be made as small as desired by choosing sufficiently 
narrow wave packets, although there is no need to improve accuracy to better than experimental resolution. 
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Table 11: Final calculation 



jfcO Real T Im T % error 

0.05 2.18499c-l -1.03160c+0 0.0982 






1 


-2 


30337e- 


1 


-4 


09325e- 


1 





0956 





2 


-1 


01512e- 


1 


-4 


61420e- 


2 





0981 





3 


-3 


46973e- 


2 


-6 


97209e- 


3 





0966 





4 


-1 


39007e- 


2 


-1 


44974e- 


3 





0997 





5 


-6 


44255e- 


3 


-3 


86459e- 


4 





0986 





6 


-3 


34091e- 


3 


-1 


24434e- 


4 





0952 





7 


-1 


88847e- 


3 


-4 


63489e- 


5 





0977 





8 


-1 


14188e- 


3 


-1 


93605e- 


5 





0965 





9 


-7 


28609e- 


4 


-8 


86653e- 


6 





0982 


1 





-4 


85708e- 


4 


-4 


37769e- 


6 





0967 


1 


1 


-3 


35731e- 


4 


-2 


30067e- 


6 





0987 


1 


2 


-2 


39235e- 


4 


-1 


27439e- 


6 





0968 


1 


3 


-1 


74947e- 


4 


-7 


38285e- 


7 





0985 


1 


4 


-1 


30818e- 


4 


-4 


44560e- 


7 





0955 


1 


5 


-9 


97346e- 


•5 


-2 


76849e- 


7 





0956 


1 


6 


-7 


73390e- 


5 


-1 


77573e- 


7 





0992 


1 


7 


-6 


08794e- 


5 


-1 


16909e- 


7 





0964 


1 


8 


-4 


85672e- 


5 


-7 


87802e- 


8 





0956 


1 


9 


-3 


92110e- 


5 


-5 


42037e- 


8 





0967 


2 





-3 


20000e- 


5 


-3 


80004e- 


8 





0966 



These tests suggest that this method can be applied to compute scattering observables at the few GeV scale. 



IX. SUMMARY AND CONCLUSION 



In this paper we introduced a formulation of relativistic quantum mechanics that uses Euclidean Green functions or 
generating functionals as dynamical input. The motivation for this approach is to formulate few-body models at the 
few GeV scale. The models are quantum mechanical which means that they are formulated in terms of linear operators 
on a model Hilbert space. They can be treated using standard quantum mechanical methods. The advantages of this 
framework over some conventional treatments of relativistic quantum mechanics is that there is a formal relation to 
Euclidean Lagrangian field theory. Specifically, the quantum theory discussed in this paper becomes the quantum 
mechanical formulation of the field theory when the model Green functions are replaced by the Green functions of the 
field theory. A second advantage of this formalism is that all calculations can be performed in the Euclidean domain, 
without using analytic continuation, however the quantities computed are ordinary Minkowski scalar products of 
normalizable vectors. This permits the use of ordinary Minkowski-space methods. 

The structure of the physical Hilbert space was discussed in section 1 in terms of a Euclidean generating functional 
for a scalar field theory. Methods to model generating functionals in terms of connected Euclidean Green functions 
were also discussed. Representations in terms of Euclidean Green functions were given in section 6. Examples 
illustrating the equivalence of the Hilbert space inner product expressed in terms of Euclidean Green functions and 
the standard Minkowski inner product were given in equations (6.4) and (6.5). 

The Poincare Lie algebra is realized by interpreting the real Euclidean group as a complex subgroup of the com- 
plex Lorentz group on the physical Hilbert space. Self-adjoint generators that satisfy the Poincare commutation 
relations are extracted by considering infinitesimal Euclidean transformations which become infinitesimal complex 
Lorentz transformations on the physical Hilbert space. Expressions for the generators were given both in terms of the 
generating functional in section 3 and directly in terms of the Green functions in section 6. 

In section four the generators introduced in sections three and six are used to construct a mass operator whose point 
eigenstates correspond to particles. The translation and rotation operators introduced in sections three and six are 
used to construct operators that create single particle states of a given sharp momentum and spin out of normalizable 
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mass eigenstates. These states necessarily transform irreducibly under the Poincare group if they non-degenerate. 

The operators that create single-particle states are used to construct two-Hilbert space injection operators in the 
Haag-Ruelle formulation of scattering. The existence of the wave operators can be checked using a generalization of 
the standard Cook's method |19j used in non-relativistic quantum mechanics. The two-Hilbert space wave operators 
are computed using time-dependent methods with the Kato-Birman invariance principle to reduce the computation 
to the evaluation matrix elements of polynomials in . The feasibility of this method for computing scattering 
observables was established using an exactly solvable non-relativistic test model. This model demonstrated that it 
is possible to perform accurate calculations over a wide range of energies using this method. Heuristic guidelines for 
choosing the "time" n, temperature (3, and the degree of the polynomial were established in the context of this model. 
Both the "time" limit and polynomial approximations were shown to be very accurate. Even though the input to 
the scattering theory involves Euclidean quantities, scattering emerges because the scattering states asymptotically 
oscillate in phase with free-particle states as the "time-paramenter" , n, gets large. The limits do not depend on the 
exponential fall-off that is used in lattice calculations. This means that calculations can be performed at relativistic 
energies without the complications that arise in alternative formulations of scattering involving Euclidean quantities 
|20j . Since the output of these calculations are wave operators, the relativistic intertwining properties provide a 
mechanism for performing finite Poincare transformations on all scattering states. In addition, this also means the 
even perturbative approximations to the Green function or generating functional will yield unitary, Poincare invariant 
approximate ^-matrices. 

Real model calculations are done using Euclidean Green functions. In principle the model Green functions should 
be Euclidean covariant, real, satisfy cluster properties and be reflection positive. Taken together these are restrictive 
conditions. While good models should satisfy these conditions, some models may fail to satisfy some of these conditions 
exactly. One may still get good approximations if small operators are responsible for the difference between the exact 
and approximate Green functions. These questions require additional study. A relevant observation is that because 
all of the calculations are done without analytic continuation, errors will be small if \Z[fi — Off] — Z exact [fi9fj]\ is 
small on a suitable domain of test functions. One does not have to worry about errors being amplified by analytic 
continuation because all observables are expressed directly in terms of Z[fi — Off] or the corresponding quantities 
in terms of Euclidean Green functions. Computational questions are also relevant. The Euclidean Green functions 
require the accurate evaluation of integrals with modest dimensions, (4(N-I)) where N is the number of particles. 
While the integrands are smooth, the integrals must be calculated accurately, particularly for scattering calculations. 
The dimensions could make a sufficiently accurate calculation of the matrix elements difficult. This also needs to be 
explored using better models. 

Another important open question that requires more investigation is the implementation of this strategy in QCD, 
where reflection positivity need not hold on states that are not color singlets. For the type of calculations discussed, 
i.e one particle states and scattering states, the two-Hilbert space injection operator has range in the space of physical 
states and one expects that reflection positivity holds on this subspace, which is all that is needed to implement the 
computational strategies discussed in this paper. One interesting observation in QCD, which is often formulated in 
terms of Euclidean Green functions, is that there are no problems in dealing with relative Euclidean times between 
physical states in the formulation of scattering theory used in this paper. This is because the scattering is actually 
formulated in Minkowski space, even though the computations are performed in using Euclidean Green functions. 

This work supported in part by the U.S. Department of Energy, under contract DE-FG02-86ER40286. 



[1] B. D. Keister and W. N. Polyzou, Adv. Nucl. Phys. 20, 225 (1991). 
[2] K. Osterwalder and R. Schrader, Commun. Math. Phys. 31, 83 (1973). 
[3] K. Osterwalder and R. Schrader, Commun. Math. Phys. 42, 281 (1975). 
[4] J. Frohlich, Helv. Phys. Acta. 47, 265 (1974). 
[5] E. P. Wigner, Annals Math. 40, 149 (1939). 

[6] J. Glimm and A. Jaffe, Quantum Physics; A functional Integral Poind of View (Springer- Verlag, 1981). 

[7] A. Klein and L. Landau, J. Functional Anal. 44, 121 (1981). 

[8] A. Klein and L. Landau, Comm. Math. Phys 87, 469 (1983). 

[9] J. Frohlich, K. Osterwalder, and E. Seiler, Annals Math. 118, 461 (1983). 
[10] F. Coester, Helv. Phys. Acta 38, 7 (1965). 
[11] R. Haag, Phys. Rev. 112, 669 (1958). 
[12] D. Ruelle, Helv. Phys. Acta. 35, 147 (1962). 

[13] M. Reed and B. Simon, Methods of Modern mathematical Physics, vol. Ill Scattering Theory (Academic Press, 1979). 
[14] H. Baumgartel and M. Wollenberg, Mathematical Scattering Theory (Spinger- Verlag, Berlin, 1983). 
[15] T. Kato, Perturbation theory for linear operators (Spinger- Verlag, Berlin, 1966). 
[16] C. Chandler and A. Gibson, Indiana Journal of Mathematics. 25, 443 (1976). 



26 



[17] H. Kroger, Physics. Reports. 210, 45 (1992). 

[18] W. Brenig and R. Haag, Fort, der Physik 7, 183 (1959). 

[19] J. M. Cook, J. Math. Phys. 36, 1011 (1957). 

[20] L. Maiani and M. Testa, Phys. Lett. B 245, 585 (1990). 

[21] V. Wessels and W. Polyzou, Few Body Syst. 35, 51 (2004), nucl-th/0312004. 

[22] E.Zeidler, Oxford User's Guide to Mathematics (Oxford University Press, Oxford, 2004). 




FIG. 1: Re ((</>|(S - 1)|^)) vs n 50[MeV] a = 325000[GeV" 2 ] 



Im <S-1> vs n 50 MeV 




FIG. 2: Im ({cj>\(S - l)\<f>)) vs n 50[MeV] a = 325000[GeV 
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Real <S-1> vs n 100 MeV 




<S-1> 100 MeV 



100 200 300 400 500 

n 

FIG. 3: Re ((01(5-1)10)) vs n 100[MeV] a = 105000[GeV~ 2 ] 
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Im <S-1> vs n 200 MeV 




<S - 1> 200 M e V 
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n 

FIG. 6: Im ((<f>\(S - l)\<j>)) vs n 200[MeV] a = 32000[GeV" 2 ] 



Real <S-1> vs n 1 GeV 
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FIG. 7: Re {(4>\{S - 1)|<£}) vs n l[GeV] a = 600[GeV" 




FIG. 8: Im ((cf>\(S - 1)|<£}) vs n l[GeV] a = 600[GeV" 



Real <S-1> vs n 2 GeV 
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FIG. 9: Re((0|(S - 1)10)) vs n 2[GeV] a = 250[GeV" 



Im <S-1> vs n 2 GeV 
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FIG. 10: Im«0|(S - 1)|0)) vs n 2[GeV] a = 250[GeV 
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Bound state contribution to Real <S-1> at l[GeV] 
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FIG. 11: Bound state contribution to Re ((^(S 1 - l)|</>}) vs n l[GeV] a = 600[GeV~ 
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FIG. 12: Bound state contribution to Im ((<j>\(S - l)\cj>}) vs n l[GeV] a = 600[GeV" 2 ] 



